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Abstract 


The  appendage  drag  of  the  surface  ship,  Joint  High  Speed  Sealift  (JHSS),  has  been  estimated  by 
various  empirical  methods  with  drag  coefficients  for  the  separate  appendage  components.  These 
empirical  methods  have  been  shown  to  be  inaccurate  when  the  interaction  between  appendage 
components  is  considered.  Recently,  several  existing  methodologies  for  the  existing  appendage  drag 
were  reviewed.  In  this  study.  Computational  Fluid  Dynamics  (CFD)  using  the  Reynolds  Averaged 
Navier-Stokes  (RANS)  equations  was  chosen  to  perform  drag  analysis  on  several  configurations  with 
different  turbulence  models.  A  breakdown  resistance  contribution  list  for  all  the  appendage  elements  is 
also  provided  in  this  study. 

Administrative  Information 

The  work  described  in  this  report  was  sponsored  by  the  JHSS  Project  Office,  NAVSEA  05D1, 
Project  Manager  Mr.  Steven  Wynn.  The  work  was  performed  at  the  Naval  Surface  Warfare  Center 
Carderoek  Division  (NSWCCD)  by  the  Computational  Hydrodynamics  Division  (Code  5700)  under  job 
order  number  07-1-2125-146-48. 


Introduction 

Several  existing  empirical  methods  have  been  used  for  appendage  drag  predictions.  However, 
there  are  several  known  problems  when  multiple  appendage  components  arc  involved  in  the  estimation. 
Previous  attempts  using  typical  drag  coefficients  and  assembly  of  individual  pieces  have  had  difficulty 
related  to  accounting  for  overlapped  areas  and  interference  among  pieces;  using  the  actual  flow  field 
velocity  as  opposed  to  just  ship  speed;  scale  effects  and  the  best  drag  coefficient  selection. 

There  are  many  commercial  numerical  tools  available,  both  empirical  and  numerical  methods,  that 
are  capable  of  providing  the  powering  prediction  for  ships  [1~4].  Recently  a  separate  study  was 
performed  at  NSWCCD  which  used  two  appendage  resistance  prediction  programs  [5,  6],  which  are 
adapted  for  PC  execution,  and  are  further  coded  with  EXCEL  spreadsheet  for  easy  parametric  studies 

[V]. 


In  the  last  decades,  multiple  Reynolds  Averaged  Navier-Stokes  (RANS)  Computational  Fluid 
Dynamics  (CFD)  tools,  both  eommereial  and  university  based,  have  been  developed  and  validated  for 
maritime  applications.  In  this  study,  one  of  the  commercial  solvers,  CFD+-I-,  was  chosen  to  perform  drag 
analysis  on  the  configurations.  CFD-H-,  developed  by  Metacomp  Technologies,  has  been  efficiently 
applied  to  numerous  applications  solving  incompressible  steady  and  unsteady  flows.  CFD+-I-  allows  for 
easy  treatment  of  complex  geometries  by  using  unification  of  structured,  unstructured  and  multi-block 
grids.  In  addition,  wall-functions  can  be  used  for  further  reducing  the  grid  size  for  a  complex  geometry. 

Four  configurations  from  fully  appended  to  the  bare  hull  have  been  included  in  this  study.  The 
CFD  predictions  will  be  used  to  compare  with  the  model  test  data  [8]  and  the  contribution  of  drag  on 
each  appendage  surface  will  be  investigated  in  this  study.  The  results  can  be  utilized  for  design  reference 
in  the  future. 
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Geometry 


Resistance  and  Propulsion  Model  5653,  representative  of  the  JHSS  baseline  shaft  and  struts  (BSS) 
hull  form,  was  selected  for  this  appendage  drag  analysis.  Four  bulb  designs.  Baseline  Bulb  (BB), 
Elliptical  Bulb  (EB)  and  a  Stem  Bow  (ST)  are  included  in  the  model  test,  but  only  the  Baseline  Bulb  is 
modeled  in  the  CFD  study.  A  photo  of  the  bare  hull  shell  of  Model  5653  is  in  Figure  1  where  the 
baseline  bulb  is  shown  at  the  middle  of  the  three  bulb  models.  The  appendages  in  the  model  test  include 
propeller  shafts,  struts,  barrels  and  rudders  which  are  shown  in  Figure  2. 


Fig.  1  JHSS  BSS  bow  design  variations 


Fig.  2  JHSS  BSS  FA,  dry-dock  photographs 


Based  on  the  model  test  Exp  5  BSS  BH  DES  fi'om  table  B18  of  reference  [8],  the  basic  sinkage 
and  trim  information  for  the  model  position  and  the  flow  speed  of  the  bare  hull  configuration  arc  listed 
in  Table  1.  Similarly,  the  position  information  of  the  fully  appended  configuration  is  shown  in  Table  2. 
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Tabic  1  JHSS:  BSS,  bow  variations,  dynamic  sinkagc  and  pitch  for  baseline  bulb  (BB)  dynamic 


conditions;  Exp5  BSS  BH  DES 


VS 

(knots) 

Sinkage 

FP  (ft) 

Sinkage 

AP  (ft) 

Piteh  Angle 
(degree) 

FP  Draft 

(ft) 

AP  Draft 

(ft) 

Full  Scale 

36 

4.27 

-0.36 

-0.28 

33.10 

28.47 

Mode  Scale 

6.16 

0.1251 

-0.0105 

-0.28 

0.970 

0.834 

Table  2  JHSS:  BSS,  bow  variations,  dynamie  sinkagc  and  pitch  for  baseline  bulb  (BB)  dynamic 

conditions;  Exp2  BSS  FA  DES 


VS 

(knots) 

Sinkage 

FP  (ft) 

Sinkage 

AP  (ft) 

Pitch  Angle 
(degree) 

FP  Draft 
(ft) 

AP  Draft 

(ft) 

Full  Scale 

36 

4.03 

-0.05 

-0.25 

32.86 

28.78 

Model  Seale 

6.16 

0.1181 

-0.001 

-0.25 

0.963 

0.843 

Some  abbreviations  for  Table  1  and  Table  2  are  defined  as; 


BSS  Baseline  shafts  and  struts 

BH  Bare  hull 

FA  Fully  Appended 

DES  Design  displacement 
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Configurations  for  CFD  Model  and  Run  Matrix 


There  are  four  eonfigurations  ineluded  in  this  appendage  drag  study  whieh  are  shown  in  Table  3. 
The  fully  appended  eonfiguration  (FA)  ineludes  ship  hull,  four  sets  of  propeller  shafts,  four  sets  of 
struts,  four  barrels  and  two  rudders.  The  next  two  eonfigurations  are  stripped  down  from  the  fully 
appended  eonfiguration  by  exeluding  two  rudders  and  then  by  exeluding  all  the  appendages.  The 
approaeh  of  stripping  down  the  model  from  the  fully  appended  configuration  provides  a  consistent  grid 
quality  for  comparison.  For  example,  when  the  rudders  were  removed  the  grid  was  modified  only 
locally,  without  altering  the  grid  in  the  global  domain.  In  this  way,  the  flow  field  in  the  global  domain 
will  remain  similar  up  to  the  rudders.  Therefore,  the  calculation  of  the  resistance  contribution  from  the 
rudders  will  be  more  accurate  without  concerns  of  grid  inconsistency. 

An  additional  configuration  was  added  to  study  the  effect  of  the  fairwater.  The  propeller  hub  has 
been  added  to  the  last  shaft  barrel  and  a  fairwater  was  added  to  the  barrel  to  form  smoother  contours  to 
the  end  the  shaft.  This  small  geometry  extension  significantly  reduces  flow  separation  at  the  end  of  the 
second  shaft  barrel,  which  should  be  carefully  counted  in  the  measurement  for  resistance  prediction. 


Table  3  Com 

puted  configurations  for  the  JHSS  baseline  bulb  (BB) 

Configurations 

Description 

FA 

Model  5653  hull  +  propeller  shafts  +  struts  +  barrels  +  rudders 

FA  -Rudder 

Model  5653  hull  +  proi>eller  shafts  +  struts  +  barrels 

FA  -  Appendage  (bare  hull) 

Model  5653  hull 

FA  +  propeller  hub  and  fairwater 

Model  5653  hull  +  propeller  shafts  +  struts  +  barrels  +  rudders  +  fairwater 

Two  turbulence  models.  Realizable  k-e  [Ref  9,10]  and  k-e-Rt  [Ref  11,12]  have  been  added  to  the 
variation  of  configurations  and  form  a  calculation  matrix  for  the  appendage  resistance  study  for  JHSS, 
which  is  shown  in  Table  4. 


Table  4  Run  matrix 


Cases  # 

Configuration 

Turbulence  model 

1 

Fully  Appended  Configuration  (FA) 

k-e-Rt 

2 

FA  +  propeller  hub  &  fairwater 

k-e-Rt 

3 

FA  -  rudder 

k-e-Rt 

4 

FA-  Appendages  (Bare  Hull) 

k-e-Rt 

5 

Fully  Appended  Configuration  (FA) 

Realizable  k-e 

6 

FA  +  propeller  hub  &  fairwater 

Realizable  k-e 

7 

FA  -  rudder 

Realizable  k-e 

8 

FA-  Appendages  (Bare  Hull) 

Realizable  k-e 
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Grid  Generation 


Unstructured  grid  topology  is  a  practical  method  to  generate  grids  for  complex  geometries.  It  is 
especially  useful  in  dealing  with  multiple  elements  in  close  proximity.  The  unstructured  topology  allows 
a  user  to  concentrate  on  building  a  surface  grid  on  each  element  without  concern  over  matching  cell 
numbers  on  the  nearby  elements.  This  grid  topology  provides  great  freedom  in  distributing  grid  cells  to 
the  most  needed  areas.  The  grid  generation  tool,  Gridgen,  which  is  developed  by  Pointwisc,  was  used  to 
generate  all  the  grids  for  this  study.  Although  Gridgen  is  capable  of  constructing  structured,  unstructured 
or  hybrid  grids,  only  an  unstructured  grid  topology  was  used  in  this  study. 

In  the  process  of  grid  generation  for  this  model,  first  of  all,  the  geometry  has  to  be  carefully 
examined  and  all  the  gaps  and  holes  in  the  geometry  should  be  repaired  before  it  is  ready  for  the  grid 
generation  tool.  The  tool  for  geometry  repair  and  preparation  is  Rhinoceros.  Second,  using  Gridgen,  the 
given  geometry  was  triangulated  to  form  surface  grid  cells  on  the  hull  and  appendages.  In  order  to 
accurately  simulate  the  flow,  grid  points  were  clustered  in  areas  where  high  velocity  gradients  are 
expected,  which  are  around  the  bow  dome  and  stem  appendages  in  this  model.  Third,  the  surface  grid 
was  extmded  outward  from  the  wall  boundaries  to  form  the  prism  layers  which  arc  needed  for  capturing 
the  detail  boundary  layer  flow  near  the  walls.  In  the  final  grid  there  are  29  prism  layers  used  in  the  near 
wall  regions.  Finally,  after  the  prism  layers  are  formed,  the  rest  of  the  domain  was  filled  with  tetrahedral 
cells.  A  sample  of  the  unstructured  grid  for  the  fiilly  appended  configuration  is  shown  in  Figure  3. 


Fig.  3  Grid  on  the  domain  boundaries  and  on  the  hull 


In  order  to  maximize  the  number  of  grid  cells  used  in  this  study  two  assumptions  have  been 
applied  to  this  grid,  which  leads  to  more  efficient  use  of  the  computer  resource.  Since  the  focus  of  this 
project  is  to  find  the  contribution  of  the  appendage  drag,  the  free  surface  calculation  can  be  eliminated 
by  applying  the  double  hull  assumption.  The  concept  of  the  double  hull  method  is  to  symmetrically 
place  a  reversed  hull  in  the  vertical  direction  to  form  an  imaginary  domain.  Therefore,  the  free  surface 
effect  is  neglected  in  the  process.  The  top  surface  indicated  in  yellow  in  Figure  3  is  the  symmetric  plane 
between  the  computational  domain  and  the  imaginary  domain.  This  symmetric  plane,  which  will  remain 
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as  a  flat  surface  during  the  calculation,  can  be  treated  as  a  rigid  boundary.  A  close  up  view  of  the  stem 
appendages  is  shown  in  Figiu'e  4. 

Furthermore,  the  ship  is  symmetrie  for  the  starboard  and  port  sueh  that  the  eomputational  domain 
ean  be  further  redueed  to  one  side  where  the  port  side  was  ehosen  in  this  study.  The  magenta  surfaec  in 
Figure  3  represents  this  symmetry  boundary  between  starboard  and  port.  The  total  number  of  grid  eells 
ean  be  signifieantly  redueed  by  using  the  eombination  of  the  symmetry  boundary  and  double  hulls 
method.  One  ean  then  eluster  grid  eells  in  the  stem  region  and  on  the  appendage  surfaees  for  better 
capturing  the  detail  of  the  flow  features. 

The  eomputational  grids  for  various  configurations  used  for  the  RANS  ealeulations  are  listed  in 
Table  5.  There  are  four  configurations  in  this  study  and  the  nomenclature  of  these  configurations  has 
been  described  in  the  previous  seetion.  In  order  to  eapture  the  detailed  flow  field  between  the  appendage 
and  the  hull  the  fully  appended  configuration  required  9. 1 5  million  cells  where  a  significant  amount  of 
grid  cells  were  used  in  the  stem  region.  Approximately  one  million  grid  cells  were  eliminated  by 
removing  the  mdder  from  the  fully  appended  configuration  in  configuration  2.  In  the  third  configuration 
for  the  bare  hull  only  2.3  million  cells  were  used.  Finally,  in  configuration  4,  fully  appended  with 
propeller  hub  and  fairwater,  the  total  cell  count  increased  slightly  from  the  fully  appended  configuration. 


Table  5  Grid  size  for  each  configuration 


Configuration 

Description 

Number  of  cells  (million) 

1 

FA  (fully  appended) 

9.15 

2 

FA  -  mdders 

8.18 

3 

FA  -  appendages 

2.30 

4 

FA  +  propeller  hub  and  fairwater 

9.30 

Fig.  4  Grid  distribution  on  the  hull  and  appendages 
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Boundary  and  Flow  Conditions 


The  boundary  conditions  for  the  computational  model  are  shown  in  Figure  5,  where  the  two 
blocks  arc  identical.  The  block  at  the  left  is  shown  solid  and  the  block  at  the  right  is  shown 
transparently,  which  allows  us  to  view  the  boundaries  at  the  back  side  of  the  block. 

The  purple  face  on  the  top  of  the  model  is  the  rigid  boundary  from  the  double  hull  method, 
which  was  defined  as  a  slip  wall.  The  blue  face  is  the  symmetric  boundary  for  geometry  symmetry.  The 
outlet  boundary,  which  can  be  identified  as  magenta,  is  located  at  the  far  right  side  in  the  transparent 
block  in  Figure  5.  The  inlet  velocity  was  specified  as  6.16  knots,  at  model  scale,  which  is  equivalent  to 
36  knots  at  full  scale.  The  red  face  is  the  pressure  outlet  and  the  orange  face  in  the  transparent  block  on 
the  right  side  of  Figure  5  indicates  the  velocity  inlet.  The  rest  of  the  boundaries  at  the  bottom  and  the 
opposite  side  of  the  blue  face  are  far  field  conditions. 

The  computational  domain  is  two  ship  lengths  from  bow  to  the  inlet;  three  ship  lengths  in  the 
downstream  direction;  one  and  one-half  ship  length  beneath  the  hull  and  three  ship  lengths  away  from 
the  hull  in  the  port  direction. 


Fig.  5  Boundary  Conditions 
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CFD  Solution  Analysis 


Predicting  the  resistance  force  and  detailed  flow  field  about  the  surface  ship  involves  solving  the 
Navicr-Stokes  equations.  The  most  practical  method  today  for  solving  the  Navier-Stokes  equations  is  the 
Reynolds  average  Navier-Stokes  (RANS)  method.  The  RANS  method  accounts  for  the  effect  of 
turbulence  by  involving  additional  equations  which  model  the  Reynolds  stress  terms  in  the  Navier- 
Stokes  equations.  Numerous  models  have  been  developed  and  applied  to  different  applications.  For 
incompressible  flow,  the  k-z  model,  a  two-equation  turbulence  model,  has  been  widely  used  and  has 
proved  effective  in  predicting  simple  wall  bounded  shear  flows.  Here  two  additional  equations,  the 
turbulent  kinetic  energy,  k,  and  the  turbulent  kinetic  energy  dissipation  rate,  e,  arc  used  to  model 
turbulence  in  the  Navier-Stokes  equations.  However,  in  some  cases  the  traditional  k-z  model  docs  not 
adequately  predict  separated  flow  with  an  adverse  pressure  gradient  in  wall-bounded  flows. 

Consequently,  two  improved  methods  for  modeling  turbulence  were  also  used  in  this  study.  One 
method  is  the  realizable  k-z  model,  developed  by  Shih  [9],  which  includes  a  new  formulation  for  the 
turbulent  viscosity  and  a  new  transport  equation  for  the  dissipation  rate.  The  realizable  k-z  model  has 
been  widely  used  in  recent  years  and  indicates  good  agreement  with  other  k-z  model  predictions  [10] 
where  appropriate.  The  second  improved  method,  the  k-e  Rt  method,  originally  developed  by  Goldberg 
[11],  carries  an  additional  equation  for  the  undamped  eddy  viscosity,  R,  which  is  independent  of  the  k 
and  £  equations.  The  k^/R  term  has  been  used  to  replace  the  dissipation  rate,  8,  as  the  dissipation  term  in 
the  k  equation.  The  leads  to  a  lowering  of  the  turbulence  level  in  recirculation  zones  enabling  improved 
prediction  of  such  regions  [12].  Originally  the  k-8  Rt  method  in  CFD+-f-  was  developed  to  resolve 
compressible  wall  bounded  flow  and  it  has  been  proven  to  do  better  than  the  standard  k-e  model  in 
solving  separated  flow  regions.  Both  the  standard  realizable  k-z  model  and  the  k-z  Rt  model  are  used  in 
this  study.  The  comparison  will  be  discussed  in  the  following  sections. 

Total  Drag  Comparison 

The  comparison  of  the  total  resistance  between  the  measured  data  and  the  CFD  calculations  is 
provided  in  Table  6.  There  are  three  configurations  in  this  comparison,  the  fully  appended  (FA),  the 
fully  appended  less  rudders,  and  the  bare  hull.  The  second  column  in  Table  6  shows  the  measured  total 
drag  in  pound  force  for  each  configuration.  The  third  and  fourth  columns  show  the  CFD  predicted  total 
resistance  with  k-z  Rt  and  realizable  k-z  turbulence  models.  The  difference  between  the  measured  data 
and  CFD  prediction  with  the  same  configuration  divided  by  the  measured  data  forms  the  deviation 
percentages  which  are  shown  in  the  brackets.  Based  on  this  comparison,  the  k-z  Rt  model  over  predicted 
the  resistance  of  the  fully  appended  configuration  by  6.7%,  it  also  over  predicted  resistance  for  the  fully 
appended  without  rudders  by  5.9%.  But  the  k-z  Rt  model  under  predicted  the  bare  hull  drag  by  2.6%.  In 
comparison,  the  realizable  k-z  model  under  predicted  all  the  configurations  by  2.6%  to  6.1%.  The  CFD 
model  in  this  study  used  the  double  hull  approach,  which  forms  no  hull  wave.  Typically  the  hull  wave 
contributes  additional  drag  to  total  ship  resistance.  Therefore  the  total  resistance  fi’om  a  CFD  prediction 
with  the  double  hull  approach  should  be  less  than  the  measured  data.  Since  the  realizable  k-z  model 
predicted  less  total  drag  than  the  measured  data,  its  prediction  is  consistent  with  the  fact  that  surface 
waves  are  not  part  of  the  simulation. 
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Table  6  Comparison  of  the  total  resistance  between  the  measured  data  and  the  RANS  calculation 


description 

Fx  (measured) 

Ibf 

Fx  (Computation)  Ibf 

[Deviation  (%)] 

k~e  Rt 

Fx  (Computation)  Ibf 

[Deviation  (%)] 

Realizable  k-e 

FA 

34.38 

36.28  [6.7%] 

33.47  [-2.6%] 

FA 

(no  rudders) 

33.41 

34.96  [5.9%] 

32.40  [-3.0%] 

FA 

(no  appendages) 

28.72 

28.23  [-2.6%] 

26.96  [-6.1%] 

The  appended  drag  contributions  from  the  rudders  or  from  the  entire  appendage  suite  can  be 
calculated  by  subtracting  the  resistances  in  the  third  and  forth  rows  from  the  second  row  in  Table  6.  The 
difference  and  deviation  percentage  are  shown  in  Table  7.  The  total  resistance  for  the  fully  appended 
configurations  remains  unchanged  from  Table  6.  The  rudder  drag  and  appendage  drag  are  divided  by  the 
resistance  from  the  fully  appended  configuration  in  each  column.  As  an  example,  the  measured  data  arc 
shown  in  the  second  column  where  the  rudder  contributes  2.82%  of  the  total  drag  and  the  appendage 
suite  contributes  16.46%  of  the  total  drag.  In  the  third  column  the  CFD  result  with  the  k-e  Rt  model 
predicted  3.64%  from  rudders  and  22.19%  from  appendages.  The  last  column  shows  the  realizable  k~e 
model  predicted  3.2%  of  the  total  resistance  from  rudders  and  19.45%  from  the  appendages.  Based  on 
the  previous  discussion,  we  found  that  the  realizable  k-e  model  shows  better  prediction  than  the  k-e  Rt 
model  in  the  current  appendage  drag  calculation. 


Tabic  7  Comparison  of  appendage  drag  with  measured  data  and  CFD  predictions 


description 

Fx  (measured)  Ibf 

[Deviation  (%)] 

Fx  (Computation)  Ibf 

[Deviation  (%)] 

k-e  Rt 

Fx  (Computation)  Ibf 

[Deviation  (%)] 

Realizable  k-e 

FA 

34.38 

36.28 

33.47 

FA 

-0.97 

-1.32 

-1.07 

(no  rudders) 

[2.82%] 

[3.64%] 

[3.20%] 

FA 

-5.66 

-8.05 

-6.51 

(no  appendages) 

[16.46%] 

[22.19%] 

[19.45%] 
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The  Effect  of  Propeller  Hub  and  Fairwater 


The  shaft  geometry  for  all  the  previous  analysis  involves  only  the  cut  off  shaft  after  the  second 
pod  which  matches  the  shaft  geometry  in  the  model  test.  However,  in  reality  the  full  scale  geometry 
should  always  include  a  propeller,  hub  and  fairwater.  Therefore,  the  CFD  prediction  could  be  more 
accurate  if  the  detailed  geometry  were  included  in  the  model.  The  geometries  of  the  shaft  with  fairwater 
and  the  shaft  without  fairwater  are  shown  in  Figure  6.  The  shaft  in  the  foreground  shows  the  propeller 
hub,  fairwater  and  a  transition  segment  between  the  shaft  p>od  and  the  propeller  hub.  The  shaft  in  the  far 
field  shows  a  cut  off  after  the  second  shaft  pod.  Note  that  the  configuration  shown  in  Figure  6  is  for 
illustration  purposes  only  because  there  is  no  configuration  with  mismatched  geometry  at  the  end  of  the 
propeller  shafts  in  the  CFD  models.  The  configuration  with  propeller  hub  and  its  fairwater  will  be 
analyzed  in  the  following  sections. 


Fig.  6  Comparison  of  surface  grid  on  the  appendages  with  and  without  propicller  extension  hub  fairwater 

The  total  appended  drag  contributions  from  these  two  geometries  can  be  broken  down  into 
rudders;  the  inboard  shaft  and  struts;  and  the  outboard  shaft  and  struts.  Each  group  in  column  1  of  Table 
8  represents  the  sum  of  the  elements  on  both  port  and  starboard.  The  fifth  row  shows  the  total  appendage 
drag  from  the  breakdown  and  the  last  row  shows  the  total  drag  of  the  fully  appended  configuration.  First 
of  all,  the  appendage  drag,  the  sum  of  resistance  of  rudders,  struts  and  shafts,  for  the  case  without 
propeller  hub  and  the  fairwater  is  6.583  pounds  which  is  about  19.66  percent  of  the  total  drag,  33.47 
pounds.  In  comparison,  the  appendage  drag  for  the  case  with  propieller  hub  and  the  fairwater  is  6.064 
pounds  which  is  about  18.39  percent  of  its  total  drag,  32.96  pounds.  By  comparing  the  total  drags  in  the 
second  and  the  third  columns,  there  is  0.519  pounds  drag  reduction  by  adding  the  hub  and  the  fairwater 
in  the  CFD  model  which  is  about  1.55%  based  on  33.47  pounds. 
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Table  8  Resistance  contribution  from  the  appendages 


Appendage  Elements 

Drag  (Ibf)  /  Contribution 

Drag  (Ibf)  /  Contribution 

Drag  (Ibf)  /  Contribution 

k-Z  (no  fairwater) 

k-Z  (with  fairwater) 

k-Z  (with  fairwater) 

CFD++  [base  on  33.47] 

CFD-H-  [base  on  32.96] 

CFD-h-  [base  on  33.47] 

Rudder 

1.467  [4.38%] 

1.510  [4.58%] 

[4.51%] 

S&S  (inboard) 

2.625  [  7.84%] 

2.356  [7.148%] 

[7.04%] 

S4&S  (outboard) 

2.491  [7.44%] 

2.198  [6.668%] 

[6.56%] 

Appendages  Drag 

6.583  [19.66%] 

6.064  [18.39%] 

[18.11%] 

Hull 

20-27  [60.65%] 

20.83  [63.21%] 

[63.78] 

Total  Drag 

33.47 

32.96 

Noticeably,  the  drag  from  the  rudders  shows  a  slight  increase  in  the  case  with  propeller  hub  and 
fairwatcr,  4.58%,  compared  to  the  case  without  hub  and  fairwater,  4.38%.  The  resistance  increase  on  the 
rudder  may  be  due  to  the  interaction  between  the  extended  fairwater  and  the  rudders. 

Similarly,  the  drag  from  both  inboard  and  outboard  shafts  and  struts  for  the  case  with  propeller 
hub  and  fairwater  are  significantly  reduced  compared  to  the  case  without  propeller  hub  and  fairwater. 
The  pressure  contours  on  solid  walls  are  shown  in  Figure  7  and  the  separation  regions,  the  iso-surface  of 
the  reverse  axial  velocity  u,  at  the  end  of  propeller  shaft  and  at  the  junctures  of  appendages  arc  indicated 
in  blue.  In  comparison,  the  separation  at  the  blunt  base  of  the  cut  off  shaft  on  the  right  of  Figure  7  is 
stronger  than  the  shaft  with  fairwater  (left  of  Figure  7),  which  allows  quicker  pressure  recovery  and 
increase  in  the  base  pressure. 


Fig.  7  Comparison  of  pressure  contours  and  the  flow  separation  on  the  appendages  with  propeller 
extension  hub  fairwater  (left)  and  without  propeller  extension  hub  fairwatcr  (right) 
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To  explain  the  phenomena,  the  pressure  eontours  on  the  rudders  and  appendages  of  these  two 
eonfigurations  are  eompared  in  Figure  8  and  the  elose  up  view  for  this  eomparison  is  shown  in  Figure  9. 
On  the  left  figure  of  Figure  9  the  pressure  at  the  end  of  the  propeller  hub  fairwater  shows  higher  pressure 
than  the  pressure  at  the  rear  end  of  the  shaft  shown  on  the  right  of  Figure  9.  The  higher  pressure  at  the 
end  of  the  fairwater  may  have  eontributed  to  the  drag  reduetion  on  the  propeller  shaft  shown  in  Table  9. 
However,  this  pressure  rise  may  also  inerease  the  pressure  at  the  leading  edge  of  the  rudder  and 
subsequently  inereases  the  drag  on  the  rudders.  As  a  result,  it  appears  that  the  reduetion  of  resistanee  on 
the  shaft  is  more  signifieant  than  the  resistance  increase  on  the  rudders  because  the  total  resistance  for 
the  geometry  with  propeller  hub  and  fairwater  is  less  than  the  total  resistance  from  the  geometry  with  cut 
off  shaft. 


Fig.  8  Comparison  of  pressure  contours  on  the  appendages  with  propeller  extension  hub  fairwater  (left) 

and  without  propeller  extension  hub  fairwater  (right) 


/ 

L _ 

% 

Fig.  9  Comparison  of  pressure  contours  on  the  appendages  with  propeller  extension  hub  fairwater  (left) 
and  without  propeller  extension  hub  fairwater  (right)  (close  up  view) 
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The  appended  drag  eontribution  ean  also  be  eompared  based  on  the  bare  hull  resistanee  in  Table 
9.  The  measured  bare  hull  resistanee  is  28.72  pounds  and  the  measured  resistanee  for  the  fully  appended 
eonfiguration  is  34.38  pounds.  The  resistanee  inerement  based  on  the  measured  bare  hull  resistanee  is 
5.66  pounds  whieh  is  about  19.71%.  In  eomparison,  the  CFD  predietion  using  the  realizable  k-e 
turbulenee  model  ean  prediet  the  total  drag  very  well.  For  the  ease  without  propeller  hub  and  fairwater, 
the  predieted  total  resistanee  for  the  fully  appended  eonfiguration  is  33.47  pounds  resulting  in  6.5 
pounds  over  the  measured  bare  hull  resistanee.  The  ratio  of  the  appendage  resistanee  with  respeet  to  the 
bare  hull  resistanee  is  22.63%  and  it  is  about  2.92%  over  predieted  than  the  measured  data  in  eolumn  2. 
Similarly,  the  CFD  predietion  using  the  realizable  k-e  turbulenee  model  shows  32.96  pounds  for  the 
fully  appended  eonfiguration.  The  predieted  appendage  resistanee  is  5.99  pounds  whieh  leads  to  a 
20.86%  eontribution  with  respeet  to  the  measured  bare  hull  resistanee.  The  deviation  from  the  measure 
appendage  resistanee  is  about  1.15%. 


Table  9  Resistanee  eontributions  from  the  appendages 


Total  Drag 

Fx 

Measured 

(ibf) 

CFD 

Realizable  k-e 

(no  fairwater) 

CFD 

Realizable  k-t 

(with  fairwater) 

Bare  Hull  only 

28.72 

FA  (fully  appended) 

34.38 

33.47 

32.96 

dFx 

5.66 

6.5 

5.99 

dFx/Hull 

percentage  based  on  exp. 

19.71% 

22.63% 

20.86% 

Deviation 

0 

2.92% 

1.15% 
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The  Detailed  Break  Down  of  Drag  Contribution  from  Each  Element 

The  drag  contributions  from  each  element  of  the  shafts  and  struts  are  listed  in  Table  10.  This 
break  down  table  is  based  on  the  fiilly  appended  result  with  propeller  hub  and  fairwater  using  the 
realizable  k-z  turbulent  model.  The  drag  in  the  list  is  for  a  single  element  only  such  that  the  total 
appendage  drag  in  each  column  of  Table  10  is  half  of  the  total  appendage  drag  in  the  third  column  in 
Table  8  accordingly.  A  list  of  all  the  appendage  elements  is  shown  in  Figure  10.  In  general,  the 
geometries  for  the  inner  shaft  set  and  the  outer  shaft  set  are  very  similar  except  the  first  shaft  pod,  where 
the  inner  shaft  pod  is  longer  than  the  outer  shaft  pod. 


Fig.  10  Appendage  elements  of  inboard  and  outboard  sets 

The  nomenclature  of  the  elements  in  the  table  are  listed  as  follows  and  the  outboard  is  arranged 
in  the  same  manner  but  started  with  ‘o’  instead  of  ‘i’: 


ipl 

inboard  shaft  pod  #1 

opl 

outboard 

shaft  pod  #  1 

ip2 

inboard  shaft  pod  #2 

op2 

outboard 

shaft  pod  #2 

ip3 

inboard  shaft  pod  #3 

op3 

outboard 

shaft  pod  #3 

ishl 

inboard  shaft  segment  #1 

oshl 

outboard 

shaft  segment  #1 

ish2 

inboard  shaft  segment  #2 

osh2 

outboard 

shaft  segment  #2 

istli 

inboard  strut  set  #1,  inner  element 

ostli 

outboard 

strut  set  #1 ,  inner  element 

istlo 

inboard  strut  set  #1  outer  element 

ostlo 

outboard 

strut  set  #1  outer  element 

ist2i 

inboard  strut  set  #2,  inner  element 

ost2i 

outboard 

strut  set  #2,  inner  element 

ist2o 

inboard  strut  set  #2,  outer  element 

ost2o 

outboard 

strut  set  #2,  outer  element 

ifw 

inboard  fairwater 

ofw 

outboard 

fairwater 
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The  drag  contribution  from  each  appendage  element  is  listed  in  Table  10,  where  the  elements 
from  the  inboard  shaft  set  are  listed  in  the  left  two  columns  and  the  elements  of  the  outboard  shaft  set  arc 
list  in  the  right  two  columns.  With  one  on  one  comparison,  several  elements  in  the  inboard  shaft  set 
behave  similarly  to  the  outboard  shaft  set,  where  the  items  highlighted  in  bold  show  very  similar 
contributions  to  the  total  appendage  drag.  For  example,  drag  contributions  of  the  inboard  item  and  the 
outboard  item  from  the  shaft  pod  #2,  shaft  #2,  strut  set  #1,  and  strut  set  #2,  are  equal  or  differ  by  less 
than  5%. 

Some  elements  show  significant  difference  between  the  inboard  and  outboard.  For  example,  for 
shaft  pod  #1,  the  inner  one  is  about  4%  larger  than  the  outboard  one,  due  to  the  difference  in  length. 
Some  differences  between  correspondent  elements  do  not  have  clear  explanations,  for  example,  the  shaft 
pod  #3,  the  shaft  segment  #1  and  the  strut  set  #2. 

The  significant  differences  in  drag  from  the  inner  fairwater  and  the  outer  fairwatcr  arc  most 
likely  due  to  the  interaction  between  elements.  The  fairwater  of  the  inboard  shaft  set  is  very  close  to  the 
rudder  and  the  fairwater  of  the  outboard  shaft  set  is  clear  from  other  elements.  The  interaction  between 
the  fairwatcr  and  rudder  increases  the  pressure  on  both  the  fairwater  and  the  leading  edge  of  the  rudder, 
shown  in  Figure  10.  Therefore,  the  resistance  of  the  inboard  fairwater  is  less  than  the  outboard  fairwater. 


Table  10  RANS  calculation  for  the  appendage  resistance  study 


Appendage  element  (inboard) 

Drag  (Ibf) 

Appendage  element  (outboard) 

Drag  (IbO 

ipl  (shaft  pod  #1) 

0.149 

opl  (shaft  pod  #1) 

0.145 

ip2  (shaft  pod  #2) 

0.158 

op2  (shaft  pod  #2) 

0.159 

ip3  (shaft  pod  #3) 

0.186 

op3  (shaft  pod  #3) 

0.173 

ishl  (shaft  #1) 

0.067 

oshl  (shaft  #1) 

0.095 

ish2  (shaft  #2) 

0.081 

osh2  (shaft  #2) 

0.080 

istli  (strut  set  #1  inner) 

0.049 

ostli  (strut  set  #1  inner) 

0.049 

istlo  (strut  set  #1  outer) 

0.038 

ostl  0  (strut  set  #1  outer) 

0.036 

ist2i  (strut  set  #2  inner) 

0.175 

ost2i  (strut  set  #2  inner) 

0.151 

ist2o  (strut  set  #2  outer) 

0.144 

ost2o  (strut  set  #2  outer) 

0,137 

ifw  (fairwater) 

0.049 

ofw  (fairwater) 

0.151 

Total  inboard  appendage  drag 

1.099 

Total  outboard  appendage  drag 

1.178 
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A  close  up  view  of  the  fairwaters  in  Figure  1 1  shows  lower  pressure  at  the  outboard  fairwater 
than  the  pressure  on  the  inboard  fairwater.  This  pressure  difference  may  explain  why  the  drag  acting  on 
the  inboard  fairwater  is  lower  than  the  drag  acting  on  the  outboard  fairwater. 


Fig.  1 1  Pressure  contours  on  the  inboard  and  outboard  fairwaters 
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Conclusions 


The  following  are  eoneluded  for  numerieal  predietion  using  CFD-H-  for  the  JHSS  appendage  resistanee. 

1.  In  general,  the  RANS  solver  CFD-H-  provides  very  reasonable  predietions  eompared  to  the  model 
test  data. 

2.  The  k-z  Rt  turbulenee  model  predicts  higher  resistanee  than  the  realizable  k-z  turbulenee  model, 
whieh  provides  predietions  eloser  to  the  measured  data.  It  should  be  noted  that  the  free  surfaee  was 
not  ineluded  in  this  analysis.  Therefore,  the  low  drag  predietions  for  the  k-z  model  are  arguably 
more  reasonable. 

3.  The  model  with  propeller  hub  and  fairwater  prediets  less  resistanee  than  the  model  with  eut  off 
shaft  pod.  The  gradual  eurvature  from  the  propeller  hub  to  the  fairwater  prevents  massive  flow 
separation  from  the  propeller  hub,  whieh  then  improves  pressure  reeovery,  and  signifieantly 
reduees  total  resistanee. 

4.  The  resistanee  reduetion  from  adding  the  propeller  hub  and  the  fairwater  should  be  eonsidered  in 
future  model  tests  and  CFD  analysis  for  better  performanee. 
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